Spin Distribution in Diffraction Pattern of Two-dimensional Electron Gas with 

Spin-orbit Coupling 



Cheng-Ju Lin and Chyh-Hong CherrQ 
Department of Physics and Center for Theoretical Sciences, 
National Taiwan University, Taipei 10617, Taiwan 
(Dated; January 19, 2013) 

Spin distribution in the diffraction pattern of two-dimensional electron gas by a split gate and a 
quantum point contact is computed in the presence of the spin-orbit coupling. After diffracted, the 
component of spin perpendicular to the two-dimensional plane can be generated up to 0.42 h. The 
non-trivial spin distribution is the consequence of a pure spin current in the transverse direction 
generated by the diffraction. The direction of the spin current can be controlled by tuning the 
chemical potential. 



I. INTRODUCTION 

Separating spin up and spin down electrons in materi- 
als is a basic but very important procedure in spintron- 
ics. It can be achieved by the generation of a spin current 
so that electrons of different spins propagate in opposite 
directions. There have been many proposals of gener- 
ating spin current. They can be roughly classified i nto 
two categories. One is by means of magnetic methodd^^ 
and the other is by pure electrical source^SHU The first 
one usually acquires loops of electric current to generate 
magnetic field to switch electron spins. Inevitably, heat 
always comes as a by-product that is notoriously unfa- 
vored, especially when many loops are squeezed in a small 
volume. Therefore, the pure electrical methods bring the 
ultimate hope for the sake of practical applications. 

All of the electrical methods to generate spin current 
share the same mechanism, namely the spin-orbital in- 
teraction (SOI). In the presence of the SOI, up spin 
and down spin are no longer degenerate at each point 
in the momentum space. The spin orientation and mo- 
mentum lock with each other in different ways for dif- 
ferent SOI systems. Due to the time-reversal symmetry, 
an electron at a particular momentum can always find 
its Kramer's partner at the momentum of the opposite 
direction. When the SOI is very strong, some systems 
form topological insulators, where there is a gap in the 
bulk bands and there are odd numbers of the Kramer's 
pairs of the edge modeP''^. 

Although SOI plays an important role in the electrical 
means to control electron spin, it is the majo r con tribu- 
tion to the spin relaxation at the same timfPSl. The 
stronger SOI is, the shorter spin relaxation time will be. 
The reason is the following. In understanding the spin- 
momentum locking, one can effectively associate a ficti- 
tious magnetic field with each point in the momentum 
space so that the spin orientation prefers to be parallel 
to the fictitious magnetic field. Suppose that an electron 
moves in a certain momentum, its spin is not necessary 
in the spin eigenstate of that momentum. Consequently, 
the electron spin precesses and the information of elec- 
tron spin is lost in the diffusive transport. Unfortunately, 
the pros and cons of the SOI are always accompanying 



with one another. Therefore, one of the most important 
questions in spintronics for real applications is to com- 
promise the advantages and the disadvantages. In this 
regard, a robust effect of spin splitting is needed for the 
spin relaxation time to be long enough for a practical 
purpose. 

In this paper, we propose a new effect of spin-splitting 
from the coherent transport by the electron diffraction. 
Electron diffraction by the quantum point contact in the 
two-dimensional electron gas (2DEG) has been observed 
by the scanning probe microscopy technique. We will 
further point out that there is a non-trivial spin distri- 
bution in the diffraction pattern if SOI is present. It can 
be measured if the spin-resolved experimental technique 
is performed, for example, Kerr rotation spectroscopy. 

The physical picture can be summarized as the fol- 
lowing. Because of the wave nature, electrons pick up 
a transverse momentum when scattered by the quantum 
point contact, which plays a role as a selector of propaga- 
tion direction. Because of the spin-momentum locking, 
electron spin passing through the quantum point con- 
tact is determined due to the selection of propagation 
direction. After diffracted, the initial electron spin will 
certainly not be in the spin eigenstate of the diffracted 
momenta. Therefore, the electron spin precesses. Elec- 
trons in different propagation directions precess in differ- 
ent ways resulting in a spin distribution in the diffraction 
pattern. Different from the spin precession in the spin- 
relaxation mechanism, the precession in our case is co- 
herent due to the coherent transport of the diffraction as 
long as there is no mixing between spin states in different 
bands at the same momentum. 

This paper serves to provide many calculation details 
in our previous Letter^. The structure of the paper is 
given by the following. The model of the theoretical pro- 
posal and the formalism will be given in the section II. 
In particular, we illustrate our method using the Rashba 
system. In section |III A[ we show the numerical results 
of the diffraction by a single slit. In section |IIIB[ we 
improve the efficiency of the spin splitting by consider- 
ing the diffraction by a grating. In section |IV[ physical 
origin of the new spin-splitting effect is provided The rel- 
evancy for being realized in experiments is discussed. In 
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the appendix, the results for the Dresselhaus system are 
given. Some formula of the special function used in the 
computation is included in the last appendix. 



II. THE MODEL 

The system we are considering is a 2DES with SOI 
diffracted by a single slit. In experiments, a single slit 
may be realized by a quantum point contact. The prob- 
lem we want to solve is the spin distribution in the diffrac- 
tion pattern. The effective Hamiltonian of a 2DEG with 
SOI is given by the following 



H 



P 

2m* 



- ^ + a{dxPy - (fyPx) + P{(fxPx - (fyPy), (1) 



where m* is the effective mass of the electron, fj, is the 
chemical potential, and cffc are the Pauli spin matri- 
ces, and a and (3 are the SOI strength corresponding 
to the Rashba and the Dresselhaus couplings, respec- 
tively. Eq. ([l]) is a 2 X 2 Hamiltonian. The energy bands 

can be obtained easily by Ep — — A* i where 
Ap = y/ (a^ -I- /32)p2 _)_ -ia(3pxPy Apparently, SOI cre- 
ates a band crossing. denotes the upper (lower) 
band. The eigenstates are v± = (e*"^,±l)^/^/2, where 



= tan-i( 



tan-i(- 



Pp^+apy 



). For the upper band, 



), and for the lower band 0_ = 



+ TT. 

The 



I3p^+ap 

(j) contains the information of the spin orientation 
orientation of the spin eigenstate is locking with the di- 
rection of momentum. Therefore, in the system with 
SOI, there is no longer spin degeneracy. Moreover, the 
spin orientations of the eigenstate v± lie in the xy plane. 
Without perturbations, the electron spin in the 2DEG 
has only xy component. Spin orientations are opposite 
for the upper band and lower band at each momentum. 

The situation we consider here is similar to the case of 
the single slit diffraction in optics. A slit of width d is 
located at y' = 0. We compute the diffraction pattern 
of the electron wave after propagating to the screen at 
distance L away in the x-direction. The diffraction pat- 
tern is the superposition of the quantum waves from the 
slit, according to the Huygens principle. The wave am- 
plitude on the screen is the quantum superposition of all 
the spherical wave emitted from the slit. Since the wave 
nature of electron is identical to photon, the diffraction 
pattern should be the same. However, it is not trivial of 
the spin distribution in the diffraction pattern when the 
SOI interaction is considered. 

We compute the quantum amplitude using the Green's 
function method pioneered by Feynman. The quantum 
amplitude of the electrons from [x' , y') at t = to (x, y) 
at time t is denoted by {x,y,t\x',y',Q). At < = 0, elec- 
trons pass through the slit. Suppose that the slit locates 
at a; = 0, and the screen is placed at a; = L away, the 



wave function on the screen is given by 

fd/2 

^P{L,y,t)^ dy'{L,y,t\0,y',0)cl){0,y\0), (2) 

Jd/2 

where (/>(a;', y', 0) is the initial wavefunction at the slit. In 
the following, we compute Eq. ^ in the pure Rashba and 
pure Dresselhaus cases. We shall use the pure Rashba 
case to illustrate our calculation procedure. The one for 
the Dresselhaus case is given in the Appdendix \K\ 

The pure Rashba case is the one for a 7^ and 13 — 0. 
Eq. M can be written hy H = KHdiag^\ where 
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The propagator in the momentum space is given by 
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The propagator in real space is the Fourier transforma- 
tion of Eq. ^ giving as the following 



{x,y,t\x' ,y' 



00 /'OO 



dpxdpy 
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pdp ■ eiP(5™'^'^f+*''»i'^p){7(i), (5) 



where x — x ~ x' , and y ~ y — y' . Eq. ([5]) is a 2 x 2 
matrix. The matrix element can be computed in terms 
of the special functions. In Appendix |Bj we include the 
useful integration formulae of the special function. Using 
Eq.flBSl) and Eq.([B2l), we obtain 



{x,y,t\x' ,y' ,0)11 = {x,y,t\x' ,y' ,0)22 
= {—) 



(6) 
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where iFi{a;b;z) is the hyp ergeo metr ic function and 



we can compute the other elements given by 
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Similarly, 
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The propagation time t is determined hy t — niL/ppx, 
where ppx is the a:-component of the momentum of the 
electron at the Fermi energy. For small angle diffraction 
as usually considered in optics, it is a good approxima- 
tion. Furthermore, we introduce the following dimension- 
less quantities for computational convenience. We define 
a = maL/h and /3 = mjSL/Ti, which are the dimension- 
less strength of the couplings, d = d/ L and y = y/L, and 
k = pFxd/2h. 

The computation can be further simplified as the fol- 
lowing. There are four terms to compute in the ma- 
trix product in Eq. ([2]). In order to make the numerical 
results more convergent, we express the hypergeometric 
functions in the integral form. Given the initial wave- 
function (t)-{0,y',0) = ((/)i,(^2) = («,-l)^/\/2d from the 
lower Rashba band, we obtain 



d/2 

-d/2 

i 



[x,y,t\x',y',0)nMO,y',0) 

nl ( ia^d 



\2mth) ^ 



- (2n)! V fc 



n=0 

xJ./d,(,-l,-<"«V.£"d.'exp(.=^, 



13 (iS^) 



(9) 



where — x"^ +y'^ —2yy' +y''^ w r'^—2yy' and — x'^+y'^ 
in the y' integration are used, and F{y) and fn{y) are 
given by 
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The other two terms can be also simplified in the same 
way. 
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Combining Eq. (|9]), Eq. (fTol), Eq. (ITT]), Eq. the wave 
— ^ ,1 — — 

function = {ipi jip-z ) on the screen can be com- 
puted as the following 



4^ [iF{y) - a{l - iy)G{y) + adH{y)\ , 



=A [-F{y) - ia{l + iy)G{y) - ladH^y)] , (13) 
where A — -k= ( „'".f^ ) es'^*. The diffraction pattern 

,/2W V 2Trith / ^ 



/2d \ 2Trith ) 

for the electrons from the upper Rashba band can be 
also obtained in the same way. Given the initial wave 
function 0+(O,y',O) = l\/2d, the wavefunction 

■0+ = (?/;J^\ t/ij^-*) on the screen is given by 

^^^^ = A {iF{y) + a{\ - iy)G{y) - adH{y)] , 

4+' = A [F{y) - lail + iy)G{y) - iadH{y)\ , (14) 



We will compute Eq. and Eq. (|14| numerically in 



the next section. I-^I^ ="pi/)i|2 _(- \il)2\ is probability dis- 
tribution of electron on the screen, which is nothing but 
the diffraction pattern. As V'l and V'2 are the spin com- 
ponents, we will show later that the spatial distribution 



of IV'il and |'!/'2| on the screen are different in Eq. (13) 



and Eq. ( 14 ) . The difference implies a nontrivial spin 



distribution in the diffraction pattern that we hope to 
investigate in this paper. 

The pure Dresselhaus case can be computed in the sim- 
ilar way. The Hamiltonian now becomes H = AHdiag^^ , 
where 
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1. (color online) The probability density of 

|2 



and \'4'2\^ in the dimensionless parameters a = 0,/3 = 0,fc = 
7.5, and d = 0.01. The diffraction pattern for spin-up electron 
is the same as for spin-down electron. Since there is spin de- 
generacy in the momentum space, the total spin of the central 
peak is zero. 



Although the energy dispersion is the same as the Rashba 
case, the orientation of the spin eigenstates differs. It 
gives rise to different spin distribution. Using the same 
computational procedure, we obtain 

/,(-) 



^P\-> = A [F{y) + P{i - y)G{v) + iMH[y)\ , 

A [-F{y) - /3(z + y)G(y) -f ipdH{y)\ , (16) 
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where "0^ ^ — {ij^i \ V'2 '') ^^'^ the wavefunction on the 
screen for the lower Dresselhaus band. 



= A [F{y) 



A [Fiy) 



P{i - y)G{y) - iPdH{y)] , 

${i + y)Giy) + ipdH{y)], (17) 



where t/)^"*"^ = (ipl^ ,ip2 ) wavefunction on the 

screen for the upper Dresselhaus band. This completes 
all of our analytical results in this paper. In the next 
section, we will solve them numerically and investigate 
the magnetic property. 



III. NUMERICAL RESULTS 

A. Single-slit diffraction 

We first plot and the spin components without 
SOI in Fig. ^. Since we work with dimensonless pa- 
rameters, this is the diffraction pattern equivalently at 
L = 1. Furthermore, we consider the wave propagating 
to the positive x direction. 1-01^ we obtained in this case 
is exactly the same as the diffraction pattern of photon. 
Dark fringes locate at the same position as the formula 
given in all textbooks of general physics. Without sur- 
prise, IV'iP and IV'iP are the same. These results imply 
that the electron spin remain lies in the xy plane af- 
ter diffraction. Actually, without SOI, upper and lower 
bands are degenerate. Given a chemical potential, both 



spin up and spin down pass through the slit. Therefore, 
the spin distribution for each component defined by 

are trivial. Namely < S^{y) for all components. 
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FIG. 2. (color online) The diffraction pattern and spin dis- 
tribution for the Rashba lower band in the dimensionless pa- 
rameters a = 0.5, k = 7.5, and d = 0.01. (a)The yellow line 
indicates the distribution of probability densities \ip\^- The 
blue line is the distribution of IV'iPi Bund the red line is the 
one of 1 1/^2 p. (b)(c)(d) are the spin distributions < S^{y) >. 

When SOI is considered, < S'^{y) > become nontriv- 
ial. In Fig. ^ and Fig. we show the results for 
the lower and upper bands in the Rashba system respec- 
tively. If the chemical potential is tuned above the band- 
crossing point, the electrons from the upper band has 
larger diffraction effect, because they have longer wave- 
length. If the chemical potential is tuned below the band- 
crossing point, the diffraction pattern comes from the 
electron of the lower band. Fig. ([2^) is the diffraction 
pattern, namely drawn in a yellow line, for the lower 
band. It is exactly the same as the one without SOI. If an 
experimental method that is not spin-resolved, for exam- 
ple scanning probe microscopy, is used, one can not dis- 
tinguish the difference between the system with or with- 
out SOI. However, the spin components I'f/'ip and |'!/'2p 
have asymmetric y dependence. Their dark fringes lo- 
cate at the same points. "The bright fringes" distribute, 
however, asymetrically. The asymmetry results in a new 
effect of spin splitting, in Fig. ^p) to Fig. ([2ji), we com- 
pute < S'^{y) > for different components. At y = 0, 
< S^{0) >=< 5*^(0) >= 0. The spin orientation on the 
screen at ?/ = is the same as the one of the initial wave- 
function at the slit. At y 7^ 0, the other two components 
start to develop in a way that is the odd function of y. 
The spin distribution for the electrons in the upper band 
is shown in Fig. (Isb) to Fig. (fstl). Since the spin compo- 



nents IV'iP a-nd |?/'2p are asymmetric in the opposite way 
to the lower band, the spin distribution distribute oppo- 
sitely. If we define the spin current as 1) =< S"^ > Vi, 
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our results imply the existence of the spin current in the 
transverse direction, namely Zy 7^ and ly 7^ 0. 




FIG. 3. (color online) The difTraction pattern and spin dis- 
tribution for the Rashba upper band in the dimensionless pa- 
rameters a = 0.5, k = 7.5, and d = 0.01. (a)The yellow line 
indicates the distribution of probability densities IV'P- The 
blue line is the distribution of jV'iPj a-nd the red hne is the 
one of i?/'2p. (b)(c)(d) are the spin distributions < S"(?/) >. 

One important feature in those results is that the 
< S^{y) > component becomes finite at y 7^ 0. One 
thing we learn in the Rashba system is that the electron 
spin lies in the xy plane. After scattered by a slit, a com- 
ponent perpendicular to the xy plane develops! Its mag- 
nitude grows as \y\ increases. Its maximum value can be 
as large as 0.42 h which is 84% of the = h/2. Another 
interesting feature in the spin distribution is that there 
exists y = 2/0 so that electron spin at y = yo and y — —yo 
are antiparallel. It is the position where < /S'''(yo) 0, 
since only < S'^(y) > is an even function in y. We fur- 
ther note that the spin distribution is independent of the 
chemical potential. It only depends on the SOI strength. 
We also compute the spin distribution for the Dressel- 
haus case. We provide the numerical results in Fig. ([5| 
and Fig. (|6|. 

The strength of the signal to detect spin distribution in 
experiments should be proportional to < 'iJj{y)\S^\ip{y) >, 
which is the product of the spin < S*' > and the wave 
amplitude. If the chemical potential is tuned below the 
band-crossing point, the electrons at the Fermi energy 
come from the lower band. In this case, the central 
maximum of the diffraction peak is spin-polarized. If 
the chemical potential is tuned above the band-crossing 
point, the electron taken part in the diffraction come 
from both upper and lower bands. Since the electron 
spin is polarized in the opposite direction between up- 
per and lower bands. In the central peak, the signal of 
spin polarization is almost zero. In some applications, 
the spin component perpendicular to the plane is useful. 
Although this new spin-splitting effect gives rise to the 
component, the signal is weak. In the first diffrac- 
tion peak, < 5*^ > reaches its maximum. However, the 
wave amplitude is much smaller than the central peak 



maximum. Therefore, new idea is needed to enhance the 
signal of 5*^. In the next subsection, we consider the 
grating device to achieve this goal. 



B. Diffraction by grating 

The optical grating is usually made of hundreds or 
thousands slits that make diffraction peaks more local- 
ized. It also increases the separation distance between 
the diffraction peaks of the photons of different frequen- 
cies. Therefore, it is one of the basic optical devices to 
separate photons of different frequencies. Here, we bor- 
row those beautiful features to design an electronic device 
to enhance the spin polarization. 

The computational procedure is the same as the one 
for the single-slit problem. The only difference is the in- 
tegration range. Here, we consider 10 slits and 20 slits 
with equal distance. The quantum amplitude on the 
screen comes from the initial wavefunction at all point 
sources in all slits. In Fig. Q, we show the results of 
the calculation for iV = 10 and A*" = 20 for the lower 
Rashba band, where N is the number of slits. The pa- 
rameters = 0.5 corresponding to the real parameter 
a'h{/3fL) = 7.62 x 10~^^eV ■ m are used and the screen 
is L = 1 jim away. We reproduce the correct diffrac- 
tion pattern for grating. The peak for A^ = 20 are much 
localized and sharper than the one for A^ = 10. 

The asymmetry of j-^i p and |V'2 P gives rise to the same 
spin distribution as the single-slit problem. It is indepen- 
dent of not only the wavelength of electron but also the 
number of slits. We will provide the explanation in the 
next section. The position of the diffraction peaks makes 
no difference from the one of photons. Different from 
the case without SOI, the spin orientation of diffraction 
peaks are different. One can tune the chemical potential 
to change the electron wavelength and the position of the 
diffraction peaks. Since < S'^(y) > is independent of the 
electron wavelength, one can tune the spin orientation of 
the diffraction peaks. In optics, it is the first diffraction 
peak used to distinguish photons of different frequency. 
Here, the first diffraction peak at positive y and negative 
y are different. At y = yo mentioned in the single-slit 
case, the spin orientation of two first diffraction peaks 
is antiparallel. When the chemical potential is tuned so 
that the first diffraction peaks locates at y = |yo|, one can 
observe opposite spin orientations on the screen, mimick- 
ing the Stern-Gerlach experiment. 



IV. DISCUSSION AND CONCLUSION 

As mentioned in the introduction, the spin-splitting 
effect proposed in the paper is due to the precession of 
electron spin subjected to the fictitious magnetic field by 
the SOI. After scattered by the slit, electron acquires a 
transverse momentum, namely the y direction. The elec- 
trons propagating in the positive y and negative y precess 
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FIG. 4. (color online) Distributions of probability density of 
IV'Pj IV'ii^j |?/'2|^for grating devices = 10 and 20 in 
the dimensionless parameters a — 0.5, k = 7.5, d = 0.01, 
and b — 0.03, where is number of slits, b — b/L, and b is 
the separation distance between nearest slits. They are the 
results of the Rashba lower band, (a) iV = 10. (b) iV = 20. 



mixing can happen by the thermal excitation. In usual 
cases, energy to cause interband mixing is around few 
meV, which corresponds to few tens of kelvin. 

A good feature of current proposal is that the diffrac- 
tion pattern is exactly the same as the one of photon even 
when the SOI is present. The diffraction nature depend 
hundred percents on the wave nature of electron. The 
new thing we discovered is the non-trivial spin distribu- 
tion in the diffraction pattern due to the spinor nature 
of electrons and the SOI. It brings huge convenience to 
control the peak of the diffraction peaks without wor- 
rying the strength of SOI. The experimental feasibility 
is already discussed in our precious Letter and will not 
repeat here. 

In conclusion, we propose a new effect of spin-splitting 
in the 2DEG in the presence of SOI. It is based on the 
electron spin precession due to the fictitious magnetic 
field generated by the SOI. Diffracted by slits, electron 
gains transverse momentum and processes in opposite 
ways in the different transverse momentum. Spin current 
ly and ly is generated for the Rashba case and 1^ and 
ly is generated for the Dresselhaus case. The presence of 
the spin current leads to nontrivial spin distribution in 
the diffraction pattern. Most importantly, the new effect 
shed lights on the potential applications in spintronics. 

We appreciate the stimulating discussion with Chi-Te 
Liang for various experimental aspects. This work was 
supported by the National Science Council of Taiwan un- 
der NSC 97-2112-M-002-027-MY3. 



in opposite directions, and hence the spin distribution 
arises. Due to the precession mechanism, the spin distri- 
bution depends on the distance of the screen from the slit. 
However, it does not depend on the electron wavelength. 
The reason is giving as the following. Changing the wave- 
length means to change the electron momentum. For 
example, if the electron has larger momentum, it takes 
shorter time to reach the screen. On the other hand, the 
magnitude of the fictitious magnetic field is proportional 
to the electron momentum. The larger the momentum is, 
the larger the fictitious magnetic field. Even though the 
electron with larger momentum reaches the screen faster, 
it also has larger rate of precession. Therefore, the fi- 
nal spin orientation becomes independent of the electron 
wavelength. 

The precession mechanism also provides the reason 
why the spin distribution does not depend on the num- 
ber of slits, namely gratings. It depends only on the 
strength of SOI and where to measure it, namely dis- 
tance of screen. Within spin coherent length, one can 
use this new effect to manipulate electron spin. Usually, 
the coherent transport can be observed more easily in the 
low temperature. In our system, the spin coherent trans- 
port can be achieved more easily if there is no interband 
mixing, that mixes electrons of opposite spins. Interband 
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FIG. 5. (color online) The diffraction pattern and spin dis- 
tribution for the Dresselhaus lower band in the dimensionless 
parameters /3 = 0.5, k = 7.5, and d — 0.01. (a)The yellow line 
indicates the distribution of probability densities [ip]^- The 
blue line is the distribution of IV'iPi sjid the red line is the 
(b)(c)(d) are the spin distributions < 5"(j/) >. 



one of |V'2p. 
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FIG. 6. (color online) The diffraction pattern and spin distri- 
bution for the Dresselhaus upper band in the dimensionless 
parameters /3 = 0.5, k = 7.5, and d — 0.01. (a)The yellow line 
indicates the distribution of probability densities Itp]^- The 
blue line is the distribution of IV'il^j ^'iid the red line is the 
one of |V'2p. (b)(c)(d) are the spin distributions < S'{y) >. 



Appendix A: Spin distribution in the Dresselhaus 
case 

The diagonalization of Eq. ([T]) with a ~ and /3 7^ 
gives H = AHdiag^^, where 



^diag ~ 
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The propagator in the momentum space is given by 



cos(fi) -ie*»sin(ft) 



''^sinift) cos(2£t) 



PP4 



U21 U22 



(A3) 



The kernels in the Dresselhaus case can be written in 
terms of hypergeometric function given by 

(x, y, t\x', y', 0)11 = (x, y, t\x' , y' , 0)22 



\2TrithJ ^ 



(2n)! V 



x,^,(n+l;l;— ) 



(A4) 



^^y. (n + 1)! / 2im/3^ t\" 



tF,(n + 2;2;— ), 



(A5) 



and 



(a;,2/,t|x',y',0)2 



f ^) e^''\ix + y) 



E 

n=0 



2nithJ \ fi / 

(n+1)! / 2ima^t^ " 



(2n+ 1)! 



i^i(n + 2;2;— ). 



(A6) 



For the electron comes from the lower band, the ini- 
tial wave function is ^_(0,y',0) = (1, -l)'^/\/2d. We 
compute the matrix product in Eq. ([2| as the following 



d/2 



x,y,t\x',y',0}nMO,y\0) 



-d/2 

1 / m 



y/2d \2TTithJ {2n)\ \ k 

1 r c^i"^ 

X i, ds(s- l)-("+i)s" / d2/'exp(s- 

27ri ./ ./_^/9 



2?ii 



1 / md \ i , , 



(A7) 



where we used an approximation = a;^ + — 2yy' -I- 
~ _ 2yy'j — + y^ in the y' integration and 
the functions 



^1 fip-^d\ , , , 



/n(2/) 



(2n)! V 



d \" ^ imr'^ ^sm{sky) 

ds J 2ht ky 



(A8) 
(A9) 



Furthermore, 



d/2 



— 1 / JTl \ i 

(L,y,^|.',y',O),20.(O,y',O) = -^(— )e. 



-d/2 V2d 
m/3\ ^ (n + 1)! A/S^JN " r(n + 1) 



E 



^^^J2n + l)! V fc ; r(n + 2) 

X — /ds(s- l)-("+l)5" + l 

27ri / 

X / dy'(-iL + y-y')exp(s— — -) 

-d/2 

1 / md \ , - 

i^(2^)^'^^*^t(-.-)G(y) + .^-i^(.)] 



(AlO) 
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where 



Giy) = E 



(2n + 1)! 



9n{v) 



d \ " imr'^ sm( sky) 

— s expfs— — — ) ; 

ds J 2ht ' ky 



s=l 



hn{v)-- 



d \ „ imr^ sky cosi sky)— sm{ sky) 



(All) 

s=l 



For we have 

l.d/2 



d/2 



{x,y,t\x,y,0)22<p2i0,y, 



J2d\2mtfi 



(A12) 



and 



(i/2 



— 1 / ?Tl \ i 

m/3\ ^ (n + 1)! J\ " r(n + 1) 

2in J 

(A13) 

where F{y), G{y), andH{y) are given by Eq. (34), Eq. 
(36), and Eq.(38), respectively. Therefore,, the electron 

wave function on the screen tp^~'^ = \ip2 ^) 
Dresselhaus lower band can be simplified as the following 

[F{y) + Pit - y)G{y) + iMH{y)] , (AM) 
=A [-F{y) - /3(z + y)G{y) + ipdHiy)] . (A15) 

Similarly, for Dresselhaus upper band 

[F{y) - P{i - y)G{y) - iPdHiv)] , (A16) 
4+) ^A [Fiy) - Pit + y)Giy) + Z;9Ji?(y)] , (A17) 



whereA = -3= (.P^w) e^^*- 



Appendix B: Special functions 

In this section, we provide some useful integral for- 
mula for the computation of the kernel in Eq. ([2|. The 
hypergeometric function iFi(a;b] z) can be expressed in 
an integral form given by 



i_Fi(a;c;z) 



ni + a-c) (c-l)!^^^^^_^^__,^,_,^,. 



r(a) 27ri 



(Bl) 



where the integral contour is around the pole s = 1 and 
c is an integer. It is involved in the integration 



n! 



u2 r°o 



2a" 



where 



^ iFi(7i+l;l;--)= / e-'^^ Jo(6y)y^"+My, 



rp 1 



h 2tt Jq 



2tc 



,^p{xcos ij>+ys 



(B2) 



(B3) 



is the Bessel function of the first kind, and f = yjx^ + y^- 
In the following, we list some formula used in the com- 
putation. 



.^^^v{x<,o^^^yAn^) ^ _^27r^J'(^) (B4) 

r n 



2tt 



)Sin 06" 



p(a: COS sin 0) _ _,-9^^ T'( — 



*27r^J(f) (B5) 
r n 



V^y'Mby)y'-+'dy^^-^^,F,in+2;2;-^J 



62 
(B6) 
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